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ABSTRACT 

Gamma Ray Burst prompt emission is believed to originate from electrons acceler- 
ated in a highly relativistic outflow. "Internal shocks" due to collisions between shells 
ejected by the central engine is a leading candidate for electron acceleration. While 
synchrotron radiation is generally invoked to interpret prompt gamma-ray emission 
within the internal shock model, synchrotron self-Compton (SSC) is also considered 
as a possible candidate of radiation mechanism. In this case, one would expect a 
synchrotron emission component at low energies, and the naked-eye GRB 080319B 
has been considered as such an example. In the view that the gamma-ray lightcurve 
of GRB 080319B is much more variable than its optical counterpart, in this paper 
we study the relative variability between the synchrotron and SSC components. We 
develop a "top-down" formalism by using observed quantities to infer physical pa- 
rameters, and subsequently to study the temporal structure of synchrotron and SSC 
components of a GRB. We complement the formalism with a "bottom- up" approach 
where the synchrotron and SSC lightcurves are calculated through a Monte-Carlo sim- 
ulations of the internal shock model. Both approaches lead to the same conclusion. 
Small variations in the synchrotron lightcurve can be only moderately amplified in the 
SSC lightcurve. The SSC model therefore cannot adequately interpret the gamma-ray 
emission properties of GRB 080319B. 

Key words: gamma ray bursts: general - radiation mechanisms: non-thermal 



1 INTRODUCTION 

Gamma Ray Burst (GRB) prompt emission lightcurves are complex with superimposed rapid short-scale variabilities. Vari- 
abilities of the order of milliseconds in the prompt phase, detected in 7-rays, were known since the discovery of the earliest 
GRBs and has led to the subsequent proposal of the "internal shock" model, where the energy in the relativistic flow is 
dissipated through multiple colli sions within the ejecta. 

In the internal shock model ( Rees fc Meszaroall994 ). the ultra-relativistic outflow from the central engine (ejecta) consists 
of a succession of shells with random lorentz factors. When a fast moving shell (with lorentz factor T / ) collides with one moving 
slowly (with T s ) ahead of it, a pair of internal sho cks develops which dissipates the kinetic energy in the flow. Each pulse in the 
burst lightcurve corresponds to one such collision ( Kobavashi et al1ll997l : lMaxham fc Zhangll2009l ). The physical parameters of 
the dissipation region, magnetic field and electron distribut ion, depend on the masses and initial lorentz factors of the colliding 
shells, and the unknown microphysics of relativistic shocks ( Daigne fc Mochkovitchlll998l ). Hence they vary erratically between 
collisions and so does the final flux. 

Almost in all cases, prompt emission has been observed only in the narrow 7-ray band until very recently. This has limited 
our understanding of the underlying emission process. Observed spectra suggest that the radiative process is non-thermal. The 
most likely candidate is synchrotron radiation. Nonetheless, sy nchrotron self-Compton (SSC) process has been also suggested 
(e.g. ( Panaitescu fc Meszarosl 200d : Kumar fc McMahon 2008)). For such models, one would expect a synchrotron component 
peaking in the lower energy band, and prompt optical emission is expected. In recent years, a few rapidly responding GRB- 
dedicated optical telescopes (e.g . RAPTOR, TORTORA, ROTSE) have become instrumental in detecting optical emission 
simultaneous to the 7-ray burst ( Vestrand et al. 20051 . 20061 ) . Most often these detections were limited to a few observations in 
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the entire duration of the burst. Nevertheless, in a few cases a temporal correlation could be est ablished between the optical 



and t he 7-ray lightcurves, indicative of their possible origin from the same dynamical process dYost et alj|2007l: Pag e et al 



20071 ). This improved spectral coverage has led to a better understanding of the prompt emission region dShen fc Zhan j 



The optical flash of GRB080319B seen in unison with the 7-ray emission was exceptionally bright (jRacusin et al 



2009). 



2008). 



Optical prompt emission was observed throughout the entire duration of the 7-ray component with remarkable time resolution. 
The onset is simultaneous in both bands and the overall shape of the lightcurves are similar, indicating that emission in the 
two bands are possibly physically related. Flux in V-band was almost four orders of magnitude higher than the extrapolation 
of the 7-ray spectrum, implying that the two lightcurves are likely to have originated from two different emission processes. 

High and low energy emission tracking each other but belonging to two different radiative processes naturally led to the 
conjecture that optical prompt emission in GRB080319B is due to synchrotron mechanism in the internal shocks and these 



photons were up-scattered to the 7-ray band by the SSC process ( Racusin et al. 20081 ; Kumar fc Panaite scu 2008)). Despite 



its advantage of interpreting the rough tracking behavior between the two bands, this model also has several difficulties. 
For example, a few seconds lag between the two lightcurves is not straightforwardly expected in this model. Several later 
calcul ations claimed that the emission radius r equired under this scenario will be much larger if internal shocks were to 
occur jKumar fc Naravanllioogl ; IZou et alJboogl ). Another drawback of this model is the energy crisis th at occurs due to the 
presence of the bright 2nd order SSC component ( Piran et al. 20091 ; Bosniak et al. 20091 ; Zou et al. 20W ). The non-detection 



of thi s second SSC bump in the prompt emission spectra as observed by Fermi LAT (e.g. lAbdo et all |2009l ); Izhang et al.l 
( 201lf )1 also places a great constraint on the synchrotron + SSC mode l. Alte rnative models to interpret the rough tracking 
optical/7-ray behavior of GRB 080319B have been proposed. iFan et al. I |2009l ) advanced the idea of a neutron loaded fireball 
where both optical and 7-ray emission are synchrotron in origin but from two diffe rent elect r on po pulations, one being the 
original electrons in the plasma while the other originates from the /3-decay process. lYu et al. I |2009l ) suggested that a pair of 
internal forward and reverse shocks could be respo nsible for the 7-ray and op tical emission respectively. Acknowledging the 
difficulty of the simplest internal shock SSC model. iKumar fc Naravanl (2009) invoked relativistic turbulence to improve the 
SSC model (cf. lLazar etafl l |200g| )). 

One inte resting observationa l feature of GRB 080319B is that its 7-ray lightcurve is much more variable than its optical 
counterpart ( Racusin et al. 20081 ) . The time resolution of optical observation is poorer than 7-rays, but even if one re-bins 
the 7-ray lightcurve to the same temporal resolution as the optical lightcurve, the 7-ray lightcurve still appears m uch more 



varia ble. This feature would give important constraints on the models. For example, the FS/RS internal shock model (|Yu et al 



2009) would predict a similar variabi lity in both the optical and 7-ray lightcurves, so it cannot interpret the above feature 



The two-zone model (|Fan et al.ll2009l ). on the other hand, is more consistent, since the optical e mission is expected to occur 
at a larger radius, w here the angular spreading time is longer. The synchrotron + SSC model ( Kumar fc Panaitescu! 20081 : 



Racusin et al. 2008j) is more difficult to access since the relative variability between the two emission components has not 



been studied in the past. 

In this paper, we study the relative variability within the framework of the synchrotron + SSC model. We approach the 
problem through two complementary methods. In the 'top-down' method, we use the observed optical lightcurve as the input 
synchrotron component, derive the fluctuations in the underlying physical parameters, and self consistently calculate the SSC 
lightcurve. In the 'bottom-up' method, we follow the standard formalism to simulate the lightcurves in the frame work of the 
internal shock model. We generate a set of basic physical parameters through Monte-Carlo simulations, calculate both the 
synchrotron and SSC lightcurves and compare the fluctuations. The aim is to compare the relative variability between the 
two lightcurves and then address whether the observational features of GRB 080319B can be interpreted. 

In section-2 and section-3 respectively, we describe our methods and results from the two approaches mentioned above. 



2 LIGHTCURVE CALCULATION IN THE 'TOP-DOWN' METHOD 

We first construct a 'top-down' method where physical parameters are expressed in terms of the observed optical luminosity. 
This approach enables us to reconstruct te mporal fluctuations of the physical parameters from the structure of the observed 
synchrotron lightcurve ( Beskin et al. 2010T ) and use it to estimate the corresponding variability that would appear in the SSC 



component. 

The synchrotron component is the input V-band lightcurve itself. We have to consistently estimate the SSC component 
that would have arisen from the optical photons and the electrons that have produced them. This would require knowledge 
of the bulk lorentz factor T of the outflow, the distance R of the emission region from the center of explosion, the co-moving 
magnetic field (£>') in the dissipation region, the electron distribution AT(7 e ), and the ratio y of SSC to synchrotron luminosity. 
The temporal structure of the observed optical lightcurve is the combined effect of the time evolutions of all these parameters. 
It is difficult but possible to disentangle each of these parameters from the optical lightcurve alone with some simplifications 
and assumptions. 
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2.1 Magnetic field in the internal shock region 

In the internal shock scenario, magnetic fields are generated by the shocks in the dissipation region (e.g., Medvedev fc Loebl 



( 19991 )). Since this process is poorly understood from theoretical considerations, empirical methods are followed where the 



magnetic energy density is assumed to be proportional to the dissipated thermal energy measured in the co-moving frame. 
If L w is the luminosity of the wind from the central engine, St is the typical variability time-scale and rj is the efficiency 
of energy dissipation, internal energy in the co-moving frame can be expressed as (r]L w St)/T, where T is the bulk Lorentz 
factor of the final shell after collision, which enters the expression through frame transformation. Assuming that the shells 
are of equal mass, F can be written as T sy /~a^, where a g = Tf/T s . If the dissipation region is at a distance R from the central 
engine, the co-moving volume can be written as AnR 2 AR', where AR' is the co- moving width of the shell. This width can 

a 2 —1 

be approximated as AR' = R/V, and the variability time scale can be expressed as St = a a < 2 . This leads to the final 

expression of the co-moving magnetic energy density u' B = B' 2 /8n = £5^77 °^ 1 r i^ e ( Zhang fc Meszaros 2002 ). 



From theoretical considerations, the bolometric luminosity is essentially related to the luminosity of the outflow (wind) 
from the central engine, L w . An amount r)L w of the original wind luminosity is dissipated as internal energy via internal 
shocks, which is carried mostly by protons. Depending on the interaction between the protons and electrons in the plasma, 
a fraction e e of this thermal energy is transferred to the random kinetic energy of the electron pool, of which a fraction k 
is radiated away. In the fast cooling regime, it is valid to assume that all the kinetic energy available to the electron pool is 
converted to radiation (k — 1). Hence the bolometric luminosity can be written as 

iboi = e e r)L w (1) 

Using equation-Q] to replace rjL w in the expression of u' B , one can express the co-moving magnetic field strength B' as 



B'^ 193./H Lboli5 J ± 1 1 (2) 

where Ri§ is R in units of 10 16 cm, T3oo is F/300 and Lboi,52 is the bolometric luminosity in units of 10 52 erg/sec. 

We require to know the bolometric luminosity Lboi to calculate B' . In the next section we describe how Lboi,52 can be 
written in terms of the observed optical specific luminosity L Vv and other physical parameters. 



2.2 From observed optical specific luminosity to bolometric luminosity 

The bolometric luminosity includes radiation emitted via both synchrotron and synchrotron self-Compton processes. Lumi- 
nosity of the first orde r IC component can be expressed as, Z/i qi = yiL syn , and that of the second order IC can be written 
as Lie, 2 = 3^i3^2L syn ( Sari et al. 19961 ; Kobavashi et al. 2007 ). where is the Compton parameter for the first order IC 



scattering defined by Lic,i/L ayn , and 3^2 is the Compton-parameter for the second order IC scattering defined by Lic,2/Lic,i- 
Compton scattering between electrons of Lorentz factor 7 and synchrotron photons of frequency v' Byn (measured in the 

co-moving frame of the relativistic ejecta) can be treated in the Thomson regime if 7 < 1. Using the characteristic 

frequency, v = 2 7rm c B'^y 2 , of a synchrotron photon emitted by an electron of Lorentz factor 7, this threshold leads to a 
limiting Lorentz factor 7kn ~ 3500(B'/1000G) -1 ^ 3 , above which Klein-Nishina corrections to the scattering cross-section 
become important. Alternatively, any SSC photon above a limiting frequency u^_ N , defined as 7KNWi e c 2 , has undergone the 
scattering process that took place in the KN regime. The Klein-Nishina limiting frequency for the first order SSC is ~ 300 GeV 
(-B'/1000G)~ 1//3 r3oo, hence the SSC scattering leading to the soft-7-ray emission can safely be assumed to be in the Thomson 
regime in the rest-frame of the electrons. As a result, one has 3^1 = 3^Th- However, following the same argument, 7kn,2 of the 
second order SSC scattering between electrons and the first order SSC photons is 117.2(_B'/ i 000G) _1,/5 and vkn,2 is ~ 18 GeV 
(-B'/1000G)~ 1//5 r3oo. Vkn,2 could fall below v^ 2 of the 2nd order IC component, which means a fraction of the scattering 
events will require KN correction. As a result, the KN correction could very well be applicable for GRB080319B, and one is 
likely to have 3^2 = [Vkn < 3^Th- Hence, the bolometric luminosity can be written as 

Lbol(*) = (1 + ^Th + J'ThJ'KNlLjnft) (3) 

Synchrotron luminosity L S yn can be estimated from the observed V-band specific luminosity Lv, once we know the 
spectral regime the optical band belongs to. A precise estimation of this requires observation of the spectral index and break 
frequencies, which we do not have for GRB080319B. Nevertheless, since we are attempting a general framework for comparing 
synchrotron and SSC variability starting from the synchrotron lightcurve, we investigate all the possible spectral regimes. 

The typical inferred magnetic field strength implies that in the prompt emission region, even the lowest energy (corre- 
sponding to a Lorentz factor j m ) electrons of the injected spectrum are undergoing heavy radiative losses. Moreover, in order 
to keep energy requirements reasonable, electrons emitting in the optical band has to be radiatively efficient. Hence the optical 
V-band is expected to be in the 'fast cooling' regime of the synchrotron spectrum. We check the consistency of this assumption 
using the value of B' obtained later in this section and find that this starting assumption is self-consistent. It is also assumed 



© RAS, MNRAS 000, [TV?? 



4 Resmi & Zhang 



that the optical band is not self-absorbed. At the end of the procedure we check the consistency of this assumption as well. 
There are only two possible spectral regimes satisfying this condition: (i) v c < i/y < u m and (ii) v c < v m < vy ■ For the former, 

the optical specific luminosity Ly can be written as L„ m ( I and for the latter it is L Vm ( 1 , where L„ m is the 

specific luminosity at v m . For fast cooling electrons, the total synchrotron luminosity L syn can be approximated as Lu m • 

r ip/2 

Hence in terms of vy, L syn can be written as Lv ^/vvv m for v c < vv < v m and Lvi/ m -^-\ for v c < v m < vv ■ 

Hence the bolometric luminosity in eq-[3]is represented in terms of the two 3^-parameters (^Th and ^kn), B' and 7 m (the 
latter two parameters entering the expression through L ayn ). 

(1 + 3^Th + y^y-a^Lyy/vyUm for v c < v v < v m 

IP/2 (4) 



(i + + yT b yKN)L v v m 



for v c < v m < vv 



From eq-[2]and eq-|4]we can see that Lboi and B' depend on each other. We essentially require B' in the rest of the formalism. 
It can be obtained by algebraically solving eq-[2] and eq-|4j and can be expressed in terms of other physical parameters. The 
final expressions of B' and Lboi are given in the appendix. The final expressions depend on the compton parameter y, which 
we will derive in the next section. 7 m also appear in the final expressions; we will be using it as the input parameter (see 
section-2.4). 



2.3 Calculation of y parameters 

In the 'top-down' method, The Compton-^ parameters enter the expression of the co-moving magnetic field through Lboi, 52- 
They will be required in the calculation of the cooling frequency as well. Before deriving these parameters, we first introduce 
the relation between the compton-parameters and the ratio e e /es which is an important equality we use throughout the 
formalism. It is used in deriving the expressions in the appendix, and also in obtaining y parameters in terms of j m . 



2.3.1 Relation with e e /es 

For the 1st order scattering, the Compton y parameter can be considered as the ratio of the SSC to synchrotron luminosity, 
which can be estimated as the ratio of the energy in seed photon field U syn to that in the magnetic field Ub- For the 2nd 
order scatteri ng, it would be the ratio between the lu minosities of the 2nd and 1st order IC components, y is equivalent to 



K Jj TjTvTW (|Sari et al. I ll99rj ; iKobavashi et al.ll2007l ) in the limit where 1st and 2nd order SSC scattering are both in the 



Thomson regime. In the fast cooling regime, where k is nearly unity, this leads to the relation, e e /es = y(l + y + y 2 )- 
However, if the 2nd order IC is in the KN regime (as is the case for our scenario in many runs), this expression is modified to 

— = 3>Th ( 1 + y?h + ^Th^KN ) . (5) 

2.3.2 For the first order scattering : 3^Th 

We first estimate 3^rh, the Compton-^ parameter in the Thomson regime, valid for the first order IC scattering. For a fast 
cooling synchrotron spectrum, if v c < vv < v m , the optical specific luminosity Lv can be expressed in terms of the peak 
luminosity L max {L v=1/c ) as, 



Lv = V Vc/vv Lmax (6) 

where L max = iV rad ^'f T (Wijers & Galama 1999). 

In terms of the total (both synchrotron and SSC) bolometric luminosity Lboi, Lv can be expressed as 

£ v _ ^bol 1_ 

(1 + yTh + ^TIi^Kn) ^JVvV m 

(here we have made use of eq-4). 

After substituting for L max , v m (= o~ ~ e -S'7m)> and v c (= „ e B'j 2 ), we obtain the relation 

= 1.95 x 10- 11 7V rad , 52 (l + ^Th + ^Th^KN)7m7cr 30 o 2 . (8) 

Since Lh °]i 2 can be substituted as 2.68 x lQ~ 5 (a g /a 2 g — l)(e e /e B )Ria^soo (using eq-2) and (1 + >rh + ^tii^kn) can be 
replaced by (l/3^Th)(e e /es) using eq-5, we arrive at, 

a, (-i -i^ — 6 ®Q -^rad.52 /^,\ 

^Th = 2.2x10 2 u2 7m 7c- (9) 

Through the same approach, the above expression can be obtained for v m < vv also. 

It needs to be mentioned that, in the standard 'bottom-up' approach, one starts from the above expression, derived by 
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Figure 1. The behaviour of 3^Th (the ^-parameter for the 1st order scattering), if the 2nd order SSC scattering is in the KN-rcgime. 
3^Th is obtained through non-linear root finding from eq-21. In the first figure, 3^xh i s plotted against 'R.rjh.E-jz for various values of e e /e.B, 
in the middle figure, 3^xh is plotted against e e /eg for a range of Ti^ASji. We can see that ^Vxh is highly sensitive to e e /es while the 
changes in IZ^AE-j^ is not very consequential. In the 3rd panel, the two y parameters, computed using eqn-21 in the appendix through 
a root-finding algorithm, are plotted for a range of x(7 m ). Here we have used an input e e /eg of 1000, 200 < f m < 230, T3oo of 1.3, 
Ri6 of 0.3 and the input optical lightcurve for Ly . For x(-~f m ) << 1, 3^KN ~ D^Thi but as i(7m) increases, the scattering enters to the 
KN-regime, the y parameter for the 2nd order scattering (in red) decreases and the y parameter for the 1st order scattering (in blue) 
increases. 



integrating the scattering cross-s ection over the electron energy spectrum assuming fast cooling electrons (jPanaitescu 
2000; lKumar fc McMahonll2008r ). substitutes for j m in terms of e e and -y c in terms of e_g, and arrives at eq-5. 



Kumar 



2.3.3 For the second order scattering : ^kn 

For the second order scattering, which is likely to be in the Klein-Nishina regime, evaluation of the 3^-parameter is more 
complex. This is because unlike in the Thomson regime, the scattering cross-section depends on the Lorentz factor of each 
electron involved in the scattering. A full numerical calculation is required to obtain an exact estimate of ^kn- In this paper, 
we adopt an analytical formalism where approximate estimates of the reduction due to Klein-Nishina effect in both scattering 
cross-section and the typical energy gain of the photons are used. Since in the SSC regime uF u peak is at v]^ for a fast cooling 
spectrum, the reduction in effective scattering cross-section can be scaled down as 



(J T 



(10) 



where x 



~trr 



7m" 



hu 



IV 



, is the normalized energy of first order SSC photon in the restframe of the relativistic electron 

300r3QQ m e c z OJ r 

with lorentz factor -y m . v 1 ^ is calculated as For x 



^> 1, one can approximate lZ a to be 



+ 



l ' 



— (log 2x^ m 

[Vkn can be written as 3^Th7^<r AZ-n, where A&r is the ratio between average energy gain for the photon in the second 
order scattering to that in the first order scattering. The average gain in the 2nd order scattering can be approximated as 

2 

7KJV 'ic"i' C ■ We divide it by the typical energy gain in the first order scattering 7^ to obtain A£n- Using the expression 
7xjv,2 = 117.2(B'/1000G) from section-2.3 and substituting for v^' 1 as 2 x 2.8 x NfiHzBy^, we can reduce AE-r, to be 
L ° Jr*™' 2 ■ Hence, when KN corrections are strong, ie., for x lm 3> 1, the final expression for ^kn is 



B ~i'Yi 



3^KN — Vl 



log 2x~, 



4 X 10 13 7KN,2 



(11) 



estimate the correction when the KN-effect is moderate (x lm is a few), hence we use the two asymptotic estimates. 

Substituting for 3^kn, we can rewrite eq-5 in terms of 3^Th, e e /es and 1Z& A£k- It is detailed in the appendix. In figure- 1 
we present the behaviour of y parameters. 

If the 2nd order scattering is in the Thomson regime, both first and second order y parameter will be the same and will 
depend only on the ratio e e /es- However, if the 2nd order scattering is affected by KN-effects, the first order ^xh will depend 
on e e /es and the term TZcAE-r.. This term signifies the extent of KN effect (For scattering in Thomson regime, it is unity. 
The more x lm is, the smaller will be lZ a AE-r. )• We estimated ^>xh for a range of e e /eB and TL a AEiz\ we find that ^xh has a 
strong dependence on the value of e e /es, and the dependence on TZaAEn is fairly weak. Only in cases of very high e e /es does 
the dependence on the 2nd term become important. In addition to that, it is easy to note from eq-8 that for a given value of 
Ee/es, a lower 3^kn (which implies stronger KN effect and larger of x(-ym)) will result in a higher 3^Th- 



2.4 Sequential steps to estimate SSC emission 

To calculate the SSC component, we require the final expressions of B' , -/V ra d, 7 m , !Vxh and ^kn (we call these class-1 
parameters), in addition to the input parameters (we will note them as class-2) a g , St, F s , p (required only for the spectral 
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regime case-ii), and e e /eB- We can see from the above sections that the five quantities in class-1 are all inter connected. The 
key is to find a way to disentangle them, start from one and to arrive at all the remaining quantities. 

A random number generator is used to determine the value of T 3 for a given pulse. Rie is computed using the relation 
of radius R = 2c5tT 2 a 2 /(a 2 — 1). Tsoo is calculated as ~ ^/ajr s /300. Across the burst, typical fractional variation in Y s , 
5Y a /Y s ~ 1/20 to 1/10. In a given simulation run a g and St are kept constant, St is kept somewhere in a range of 0.5 - 2.0, 
and a g is fixed around 2 — 5. In a given run, their values are chosen such tha t the variation in for the burst could be up 
to a factor of 3 to 4, to reproduce the range of radii in which shocks occur (jPaigne fc Mochkovitcblboooh . Between various 
realizations of the simulation, T s varies from 50 to 500. The resultant range of T3oo through multiple runs is ~ 1 to 5. Multiple 
simulations resulted in a large range of Ri6, from 0.01 to ~ 1. 

Since 7,„ is roughly constant for a given pulse during shock crossing, we use 7 m to be the input parameter, and randomly 
assign its value for a given pulse. The range of 7m varied from ~ 50 to ~ 500 through the runs. Whenever the optical frequency 
fell below v a , we moved to another set of parameters. Depending on the other parameters, e e /e_B was varied such that ^Th is 
around 10. 

Eq-[2T] connects ^Th with B' and j m . On the other hand, B' , as we can see from its detailed expression given in appendix, 
depends on 7 m and 3^Th- Other than that, only the quantities in class-1 enter in these two expressions. Hence, we can use the 
expressions of B' in the appendix and write eq I21l in terms of 7 m alone. However, we are not presenting this long expression 
(in fact one each for the two spectral regimes, vy < v m and Vy > v m ). Below are the sequential steps we follow in estimating 
the SSC emission. 

step-1: 3^Th is obtained from eq-21 for a given 7 m and class-2 parameters using non-linear root finding algorithms. 
step-2: Now we know both 7 m and 3^Th, so B' can be calculated. 

step-3: Knowing B' , 3^Th and 7m, we can now calculate 3^kn from the expression in the appendix. 

In the code, the possibility of the 2nd order scattering to be in the Thomson regime is checked by monitoring the value 
of x-y m once B' is estimated, and in that case the quantities are re-estimated. 

step-4: Once we know ^Th, !Vkn and B' we can calculate j c . If j c falls below unity, we set it to unity. 

Since we know u m and u c , we can now do a self-consistency check about the synchrotron spectral regime, and see whether 
our initial assumption of fast cooling and the location of the optical frequency is valid or not. If the consistency is violated, 
we redo the calculations with a different set of parameters (especially 7™ values). 

step-5: Finally, we use eq-[9]to estimate the value of -/V ra d with 7m, 7c and ^Th- 

We thus have all the quantities in class-1 known, and can compute the SSC lightcurve. The characteristic frequencies are 
estimated as uj: c — 2'y 2 v c and v 1 ^ = 2-y^ n u rn , and the spectral normalization at v — v' c c is estimated as I/„, m ax^Th/(7m7c). 
We follow lGupta fc Zhanj |2007h in calculating the SSC spectrum. We assu me the flux to decay following the curvature effect , 
by (t/f K )~ 2+/3 (/3 being the spectral index), after the peak (t p ) of each pulse ( Fenimore et al. 1996; K umar fc Panaitesculboool ; 



Zhang et allbood ; iLiang et al]|2006l ). 



2.5 Measure of Variability 

An important point to note is that for GRB0803 19B the time resoluti on of the BAT lightcurve is almost ~ 20 times higher than 
that of the V-band lightcurve from TORTORA ( Racusin et al. 20081 ) . Many individual pulses and fine features may have been 



lost in the coarse resolution of the optical lightcurve. Therefore a direct comparison between our synthesized SSC lightcurve 
and the BAT lightcurve is not entirely appropriate. We can only make a qualitative study on the extent of amplification 
possible by the SSC process. Our aim is to see how small-scale fluctuations of the synchrotron lightcurve would appear in the 
SSC component, and how variations in the electron distribution can control their appearance. 

The template for our analysis is the optical V-band lightcurve of GRB080319B. However, since our aim is to see whether 
the high variability of its 7-ray lightcurve is due to enhancement of small scale synchrotron fluctuations by Compton up- 
scattering, we synthesize miniscule fluctuations by adding random Gaussian fluctuations at mid-points after interpolating the 
lightcurve. We synthesize multiple lightcurves using Gaussian distributions with different widths. We first test our method in 
the original optical lightcurve from TORTORA and later use the new synthesized lightcurves. In figure-2, we present one of 
the template lightcurves with synthesized fluctuations along with the original lightcurve from TORTORA. There are a total 
of eight pulses in the lightcurve, of which three are too miniscule that they would not have been discerned as independent 
pulses typically. Observations exist during the tail only for the 3rd, 5th, 7th and 8th pulses. To obtain the beginning of the 
pulse to, we fit the decay-part of these pulses with the temporal profile expected from the curvature-effect, which requires 
the tail to fall as (t — to)~ ( - 2+/3 \ For pulse-8, which has good sampling over the tail, we could constrain the values of both 
to (18.0 ± 14.0 and the index (—2 — j3 = —2.5 ± 1.1). For the other three pulses, we have assumed /? to be around 0.5 and 
obtained the value of to- For those pulses for which we could not run a fitting routine, we used the same value of to as its 
nearest pulse. It is to be noted that the optical spectral index can be much steeper, especially because of the large spectral 
index observed in 7-rays. 

7m is the parameter that we tune to control the pulse amplification in SSC. For the above synchrotron lightcurve, there 
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Figure 2. The input synchrotron lightcurve. Left: The original data from Beskin et al. is shown in open squares connected by gray line. 
In black circles are the synthesized random fluctuations. Right : The eight pulses as we define them presented in different colors. The 
tails arc computed using the curvature effect and are shown in dash-dotted line. See text for details. 



will be a set of eight 7 m values as input. The spectral regimes where the optical and the gamma-ray bands fall play a crucial 
role in the nature of the output gamma-ray lightcurve. Given the same distribution of T3oo and Rie, we see that if i/y is 
between v c and v m and if v 1 is between i/' c and i/^p, the two lightcurves become nearly identical (even with large variation 
within the input j m set) because of the same functional dependence both synchrotron and SSC flux has on j m . This is true 
for the case where u m < vy and v]^ < as well. The best contrast is possible if vy < u m , but v., > v]^ or u~j < ul c . This is 
because of the additional dependence the SSC flux has on 7 m . 

It is easy to see, using the expression for characteristic frequencies, that vy will be below v m if B'7m > 2 x 10 11 . Similarly, 
v 1 ^ will be below 650 keV if B'y^ < 5 x 10 6 . Since in our formalism B' depends on j m and other input parameters (class-2 
parameters), this basically boils down to a range of j m for given class-2 parameters and the input synchrotron luminosity Ly. 
But since the y m dependence of B' comes through the non-linear expression in terms of 3^Th, we cannot give an analytical 
expression describing the range. 

Amplification or quenching of the SSC pulse will depend on the input j m of that pulse. An analytical expression for 
^p SC ~ 7m relation is difficult to obtain because of the non-linear nature of the algorithm. We observe the variation of Lp SC 
for a large range of 7 m values and see that the pulse peak L p is roughly proportional to 7m 3 . Hence, the ratio between two 
adjacent peaks Lp/Z/p +1 oc ('Jm/'Ym' 1 ) 3 ' 3 if Tsoo and Rie do not vary much between pulses. Ideally, any contrast between 
adjoining pulses can be produced by changing 7 m , but the pulse profile will appear unnatural with sharp dips and rises for 
stark contrasts. The variability of SSC lightcurve depends on the distribution of 7 m through pulses. In figure-3 we present a 
sample of 7-ray lightcurves we obtained along with the input optical synchrotron lightcurves. In the first panel, we have also 
given the output synchrotron lightcurve along with the input lightcurve (the only difference is the decay determined by the 
curvature-effect in case of the synthesized lightcurves) for comparison. 

In order to have a quantitative measure of the 'variability' of the lightcurves, we estimate the standard deviation of the 
burst from an average profile that best imitates the burst profile. We construct a 'trapezoidal' function with a rising part, a 
plateau and a tail for this purpose. The function is determined by five parameters: normalization /0 and nodal points a, b, c 
and d. The rising part of the function is given as f(t) = f0(t — a)/(b — a), the plateau is the constant f(t) = fo, and the 
tail is f(t) — f0(d — x)/(d — c). After each run, we fitted the entire output SSC lightcurve with the trapezoidal function by 
varying the five parameters mentioned above, and estimated the standard deviation from the best fit trapezoid. We compare 
this value with the standard deviation obtained for the optical lightcurve, which is obtained by following the same method. 

For a given distribution of Taoo and R\§ across the pulses, the best SSC amplification occurs when v c < vy < v m 
and 1*1° < v]^ < u~/. The spectral regime u-y < i^F < v]^ can also produce similar amplification, but for ul c to be below 
~ 650 keV one requires very small (~ 1 - 10 G) co-movin g magnetic field . Only alternative scenarios involving quick decay 
of magnetic fields in the shock downstream (for example, iPe'er fc Zhand l|2006l )) can achieve such low magnetic fields for 
standard input parameters. Moderate amplification can be obtained also if v c < v m < uy and ul° < v~, < v]£ . But this 
spectral combination also requires the co-moving magnetic field to be too low. We could not find any condition where large 
variability amplification is possible. Hence we conclude that variabilities in the synchrotron lightcurve can be moderately 
amplified in the SSC lightcurve. 
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Figure 3. In the first panel, we have plotted the input synthesized optical lightcurve (black) and the simulated lightcurve (red) from 
one of the runs. The rest of the panels (2 to 6), with normalized synchrotron luminosity in red and SSC in blue, showcase a variety 
of SSC lightcurves from our simulations. 7 m for the individual pulses and the values of e e /es & n d St for each runs are included 
in the panels. Computed values of r 3 oo and iJig for each runs are listed below. The numbers correspond to individual pulses. For 
panel 2 & 3 : r 300 = {0.759, 0.85, 0.787, 0.8, 0.839, 0.717, 0.984, 0.86}, R 16 = {0.124,0.144,0.134,0.119,0.122,0.121,0.18,0.128} , for 
panel 4 : r 300 = {0.83, 0.778, 0.787, 0.801, 0.839, 0.716, 0.984, 0.859}, R 16 = {9.9,8,8.9,7.9,8.2,8.1,12.0,8.6} x 10~ 2 , for panel 5 & 6 : 
r 3 o = {1.4, 1.42, 1.4, 1.47, 1.51, 1.37, 1.48, 1.5}, Ri 6 = {0.42, 0.4, 0.42, 0.4, 0.4, 0.44, 0.39}. The only difference between panel 5 and panel 
6 is the location of Uy in the SSC spectrum. In 5 v-y > while in 6 it is the opposite. In all panels, v c < vy < v m . Both vy and v 1 
being in the same region of the spectrum diminishes the variability. 



3 THE 'BOTTOM-UP' APPROACH AND SIMULATED LIGHTCURVES 

To understand the variability better and also to have a complete picture, we approach the problem from the opposite end. We 
simulate the burst lightcurves from the 'bottom-up' method where the final flux is built up from the input physical parameters 
of the emitting plasma. We assume the internal shock framework where the burst lightcurve is the sum of several independent 
pulses produced from random collisions between shells. We use a toy model where the temporal structure of a single pulse is 
determined by the evolution of the wind luminosity alone. H owever, more detailed modelli ng of the internal shock dynamics, 
including hydrodynamic simulations have been done in past ( Daigne fc Mochkovitch 2000j ) . 



We consider an isotropic wind luminosity L w (t) which varies between to and t p as a shallow (as t s with 5 between 0.5—1.0) 
function of the observed time t for a given pulse (ie., for one collision), to is the time of collision and t p is the shock crossing 
time (also the pulse peak). The total number of radiating electrons N e is jl L w (t')dt' / (Tm p c 2 ). In normalized units, 

iVtot.ss = t^tttTy^ — (12) 
450(f) + l)r 30 o 

The co-moving magnetic field (B 1 ) is calculated by assuming that a fraction eg of the shock created thermal energy 
(8pm p c 2 ) will be carried by the magnetic field, where Bp is the internal lorentz factor of the shocked shell. The downstream 
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magnetic field density u B i = ^-Opm p c 2 , where N p , the number of protons ejected ~ N e , and the co-moving volume V' = 
47rr 2 Ar, with the thickness Ar approximated as cir s h. This leads to 



The minimum lorentz factor 7 m of the shock accelerated electron distribution is assumed to be 

7m = e e (m p /m e )9 P (14) 

where e e is the fractional energy carried by these electrons. 

Total number of radiating electrons, iV ra d, is different from Atot if the plasma is in 'severe fast cooling', where time-scale 
(t'-y = i) for electrons to cool down and loose all their kinetic energy is less than the dynamical time scale. In that case, a non- 
negligible fraction of electrons would pile-up at 7 of unity and would not be available in the relativistic pool to radiate via the 
non-thermal processes. This fraction keeps on increasing as the source ages. An exact estimate requires solving the continuity 
equation involving electron injection and radiative losses, which is beyond the scope of this paper. When i 7= i <C t' dyn ('severe 

fast cooling'), the fraction of electrons remaining in the power-law can be assumed to be J =1 ■ Hence, the total number of 

'dyn 

radiating electrons, N rad = 0PLAtot where, 

PL = Min(l, (15) 

\ d y n / 

tty = i is the time scale measure in the co-moving frame for an electron to cool down to 7 ~ 1, which can be roughly expressed 

(the derivation is given in the appendix). It is relevant now to derive the expression for the 



1 



cooling break in the electron spectrum j c . 7 C — > 1 in the 'severe fast cooling' regime. During the normal fast cooling, it can 
be expressed as *7 =1 , which is same as 0pl- Hence we have the following : 



7c = 



</>PL if £7=1 3> t'dyn 

1 if t 1= \ <C idyii 



JV rad _ J 1 if £7=1 2* t'dyn 

Ntot [4>PL if *7=1 <SC t'dyn 

The Compton y parameter for the first order SSC scattering is calculated as 10 _6 A? ra d l 527m7c/-Ri6 (eq-© since A^ad and 
B' are known, y parameter for the 2nd order SSC scattering is estimated by eq-llll 

We calculate the lightcurves in both slow and fast cooling cases, unlike in the previous approach. The synchrotron and 
SSC spectral breaks and peak flux are now calculated following the standard procedure. The synchrotron and SSC spectrum 
are calculated as piece- wise powerlaws, as described in section-2. After the pulse peak (shock crossing time), we calculate the 
flux decay using the curvature effect. 

We ran a Monte-carlo simulation where the luminosity normalization, temporal index 8 and 9p are considered as random 
variables with a uniform distribution. The luminosity normalization and dp range between the typical values of 1.0— 15.0 and 
3.0 — 10.0 respectively. Either from observations or theory, there is not a clear idea about the values S. A time independent 
wind luminosity can not produce a rising, fast-cooling synchrotron pulse, hence we chose a shallow range of 0.5— 1.0 for S. We 
also consider the shock crossing time (or pulse peak) to be distributed randomly between 0.5 — 2.0. Each pulse thus generated 
are later shifted by a random to value. The final burst profile is the sum of all these individual pulses (see figure- l3~T]l . Like in 
the 'top-down' formalism, we scan the typical ranges of the parameters. R\q is scan from 0.01 to 1.0 and I^oo from 0.3 to 3. In 
a single simulation, the radius and V are changed by a small factor. A range of 0.01 — 0.1 in e e and a range of 10" 4 - 5 x 10" 4 
is scanned in ep- In a single run, these parameters are kept fixed. Value of p is fixed at 2.2. Even though the typical temporal 
index for the rising phase of an individual pulse is 35/4, if the observed frequency is between v c and v m much sharper observed 
rise profiles are possible due to the cumulative effect of multiple pulses. 



3.1 Measure of variability 

We estimate the 'variability' following the same method as in the 'top-down' scenario. We obtain the standard deviation of 
both synchrotron and SSC lightcurves from the best-fit trapezoid. We find that the SSC lightcurves are at best only slightly 
more variable than the synchrotron lightcurve. With some set of parameters, the SSC variability can even be less than the 
synchrotron one. This is because the overall flux variation in the SSC is much lesser than its synchrotron counter part. The 
number of extra 'spikes', even if present in the SSC, do not contribute much in the estimated variability. 

The luminosity normalization and dp have more influence in the amplification of the pulses compared to 5. We made 
runs by keeping S as a constant and also by allowing it to vary. Low variation in both L w and dp and constant S can result 
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Figure 4. Typical synchrotron (left) and SSC (right) burst lightcurves from the simulations. The individual pulses arc also shown. The 
typical time resolution we have is ~ 0.01 sec in the observed frame, hence that is limit of the sharpest features from the simulation. In 
this simulation, we have used e e = 0.01, eg = 10~ 4 and p = 2.2. I^oo was varied from 1.0 to 2.5, Rig from 0.08 to 0.1, 9 p from 3.3 to 
5.5, and S from 0.5 to 0.9. L m ^2 was kept at 15(t/lsec)' 5 . 

in nearly identical pulse shapes between synchrotron and SSC components. We did not find any marked difference between 
results from slow and fast cooling cases. 



4 SUMMARY AND DISCUSSION 

In this paper, we have done a comparative study of lightcurve variability of the synchrotron and SSC components of GRB 
prompt emission. Starting from a template synchrotron lightcurve with small scale fluctuations, we trace back the magnetic 
field of the emitting region and self-consistently calculate the SSC lightcurve using assumptions for the electron distribution. 
We investigate how the small scale fluctuations in the SSC lightcurve are related to the initial variabilities present in the 
synchrotron lightcurve. A multitude of temporal structures can be obtained for the SSC lightcurve depending on the parameters 
of the electron distribution function. Degree of modification of the input fluctuations are different for different spectral regimes 
within the SSC spectrum, due to varying sensitivity to y m . Miniscule changes in the electron distribution, which appear as 
indiscernible pulses in the synchrotron lightcurve, can get moderately amplified in the SSC lightcurve. Hence, in general, the 
SSC lightcurve can be more variable compared to the synchrotron lightcurve, however, not to a large extent. We complement 
the formalism with a 'bottom-up' approach where the synchrotron and SSC lightcurves are calculated through a Monte-Carlo 
simulation of internal shock model. Only moderate amplification to the variabilities could be obtained in this approach as 
well. 

We apply our method to the "naked-eye" GRB 080319B. This burst has been interpreted within the framework of the 
synchrotron + SSC model ( Kumar &: PanaitescvjEoOo] : Racusin et al Several difficulties have been raised for this model, 

including the time lag between the 7-ray and optical lightcurves, as well as the energy budget crisis from the second-order 
SSC. Here we apply another criterion, the relative variability between the two components, to investigate the validity of the 
model. We create a synchrotron template based on the observed optical data of GRB 080319B (with minor modifications to 
test how small variabilities are amplified), and calculate the expected 7-ray variabilities in various spectral regimes. We found 
that the model lightcurves are all smoother than the observed one. We then conclude that the optical/7-ray lightcurves are 
difficult to account for within the simplest synchrotron/SSC model. 

This conclusion is reached based on our analytical formalism presented in this paper. For an analytical treatment to be 
possible, we have to adopt some approximations and/or make some assumptions. We conclude by listing them and commenting 
on their validity. 

• The total bolometric luminosity is a fraction of the power dissipated by the internal shocks, which in turn is a fraction 
of the wind luminosity from the central engine. Such a treatment has been adopted in all internal shock model calculations 
in the past. 

• The ratio e e /es, indicating the fractional energy in electrons and in the magnetic field respectively, is assumed to be a 
constant. Numerical simulations have started to derive e e and es from the first principle, but none have revealed how these 
parameters depend on shock parameters. In principle, both values may evolve with shock parameters. However, since little 
understanding is achieved regarding such an evolution, we take the simplest assumption that the ratio e e /es is constant. We 



© RAS, MNRAS 000,[T]-?? 



Prompt Emission Variability 11 



note that in afterglow modeling, the assumption of a constant e e and es throughout the deceleration phase seems to fit the 
observed data well. 

• The bulk Lorentz factor T and the emission radius R are allowed to randomly vary from pulse to pulse. By doing so, we 
have assumed that multiple collisions happen at different radii between shells of varying lorentz factor. This variation did not 
significantly affect the relative variability contrast between the synchrotron and SSC lightcurves. 

• In both 'top-down' and 'bottom-up' models, we have not considered the evolution of R and V in one collision. 

• In the 'top-down' method, the underlying electron distribution is approximated to be in the fast-cooling regime, i.e. 
even the cooling time of the electrons with the lowest injection energy is shorter than the dynamical time scale. Such an 
approximation is found self-consistent given the typical value of B' we derive. In the 'bottom-up' method, both slow and fast 
cooling conditions are tested. 

• The power-law index of the injected spectrum is assumed to be 2.2, a typical value for relativistic shocks. Our results do 
not significantly depend on this value. Changing it to a different value would not affect our conclusion. 

• The synchrotron spectrum i s approximated as a m ulti-segment broken power-law representing different spectral regimes. 



The analytical treatment of SSC l|Gupta fc Zhang||2007l ) is adopted. In the expression for the peak luminosity of the synchrotron 



spectrum, we have neglected a correction facto r that takes care of the contribution from the power-law electron distribution 



( Wiiers fc Galama 19991 ; Gupta fc Zhang 2007 ). The SSC peak luminosity is assumed to be 3^Th times the synchrotron peak. 
All these are standard analytical treatments of synchrotron and SSC emission spectra. More realistic treatments would not 
affect the variability contrast in the synchrotron and SSC lightcurves. 

• For the second order scattering in the Klein-Nishina regime, the ^-parameter is approximately estimated by scaling down 
its value for the first order scattering in the Thomson regime by a factor invoking the ratio between the KN cross section and 
the Thomson cross section at the electron injection energy. 

• The assumption of a fast cooling spectrum leads to the definition of the bolometric synchrotron luminosity 

An accurate treatment would introduce a correction factor of the order of unity, but would not affect our conclusion. 

We thank the anonymous referees for constructive comments that improved the quality of the paper. This work is supported by 
NASA NNX09AO94G, NNX10AD48G and NSF AST-0908362. LR acknowledges support from the French Agence Nationale 
de la Recherche via contract ANR-JC05-44822. 
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APPENDIX : DETAILED EXPRESSIONS 

For v c < vv < u m , 

/ \ 2/3 / 2 i \ 1/3 / .. \ 4/3 

— — "'{A) (^) S(S^) 

/ \ 1/3 / 2 -i \ 2/3 , N 2/3 

B '^^ v '{w^) (^) f^feisk) <"> 

For u c < vv < v m , 

22/21 



and 



and 



L b ol,52^0.08Lv,36 ^ 5 . 45xl01 4j {—) r B ^^{ — ) 

.10/21 / w V 1/21 fal-lY 3 1 /500\ 2 / 21 



(18) 



After substituting for 7kn,2 as 500B' (see section-2 about how we obtained this expression) we can rewrite ea - 1111 as 

,/ \ -6/5 



(20) 



where x lm defined in section — in the same normalized units is ~ 5jjo (^) 5 - 
Hence the equality eq-[5]can be written as 



^xh(i + ^ Th + A(rf m , B')y Th 2 ) = ^ (21) 



where the term 0.3 ( (ioo) sir — o 6 2a; 7m + §) which essentially is the product lZ a x -£jz is written as 

A{^ m , B'). In case if the 2nd order IC scattering is in the KN-regime (which can happen for high values of y m ), A will 
be unity as both lZ a and AS-jz will be reduced to unity. 



APPENDIX : CALCULATION OF ELECTRONS IN THE 'COOLING PILE' 

If the radiative cooling is severe, a fraction of electrons cool down, loose all their kinetic energy within the dynamical timescale 
and pile up around 7 = 1. We estimate the remaining fraction in the powerlaw as -f 1 ^- x Atot, where t T= i is the time scale 

'age 

for an electron to cool down to 7 = 1. t 7 =i can be estimated from the cooling rate as follows: 

The radiative cooling is due to both synchrotron and IC processes. The total energy loss rate ^ can be written as 
(4t) S yn + (l«)ic' The second term include two components, i.e. energy loss due to 1 st and 2 nd order IC, which can be 
approximated as y (^) s n and }^kn (^) s 1 respectively. The synchrotron loss itself can be approximated as (^) g n = 
As a result, the total energy loss of an electron can be written as 



- 1.3 x 10- 9 B' 2 7 2 (1 + y Th + ^Th^KN) , (22) 



dt 

where the constant 1.3 x 10 -9 is from ,. aT 

D7rm e 



Ideally ^kn is a function of the electron lorentz factor 7 and magnetic field is a function of time. In our treatment, we 
are focused on an average estimate. We use the expression derived for and ^kn in Sect. [2~3l and hence take out the 7 
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dependence of the coefficients in eg - 1221 We use the magnetic field at the peak of the pulse as a representative value. We can 
then solve the differential equation under the assumption that the magnetic field is a constant over time, which is not a bad 
assumption as long as the cooling time scale is much shorter than the dynamical time scale. 

/i\ TO (O'Z') 

~ 1 + 1 .3 x 10- 9 (1 + y Th + y T h^KN)7oB ,2 t 
Under heavy radiative loss, the term unity in the denominator can be neglected, so that one has y(t) ~ —z — ,„ Q , 1 , — ^ — t^tj 

Substituting for 7(4) = 1, one can then obtain the approximate t 7= i = x 1CP 9 (1 + 3^rh + ^tii^kn) £>' 2 J 
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